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Abstract. We critically examine a scenario for the enrichment of the interstellar medium (ISM) in which supernova ejecta 
follow a long (10* yr) journey before falling back onto the galactic disk in the form of metal-rich "droplets". These droplets 
do not become fully mixed with the interstellar medium until they become photoionized in Hii regions. We investigate the 
hypothesis that the photoionization of these highly metallic droplets can explain the observed "abundance discrepancy factors" 
(ADFs), which are found when comparing abundances derived from recombination lines and from collisionally excited lines, 
both in Galactic and extragalactic Hii regions. We derive bounds of lO'^^-lO'"" cm on the droplet sizes inside Hii regions in 
order that (1) they should not have already been detected by direct imaging of nearby nebulae, and (2) they should not be 
too swiftly destroyed by diffusion in the ionized gas. From photoionization modelling we find that, if this inhomogeneous 
enrichment scenario holds, then the recombination lines strongly overestimate the metallicities of the fully mixed H ii regions. 
The abundances derived from collisionally excited lines also suffer some bias, although to a much lesser extent. In the absence 
of any recipe for correcting these biases, we recommend the discarding of all objects showing large ADFs from studies of 
galactic chemical evolution. These biases must also be kept in mind when comparing the galactic abundance gradients for 
elements derived from recombination lines with those derived from collisionally excited lines. Finally, we propose a set of 
observations that could be undertaken to test our scenario and improve our understanding of element mixing in the ISM. 
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1. Introduction 

The detailed process of enrichment of the interstellar medium 
(ISM) by products arising from supernovae explosions and stel- 
lar winds is far from being fully understood (see review by 
Scalo & Elmegreen 2004). Many mechanisms are likely to be 
at work on different scales, both in time and in space (see, e.g., 
Bateman & Larson 1993; Roy & Kunth 1995; Tenorio-Tagle 
1996; de Avillez & Mac Low 2002). 

The present-day chemical composition of the ISM is de- 
rived from the analysis of the chemical abundances of H ii re- 
gions. In such analysis, it is always assumed that Hii regions 
are chemically homogeneous. So far, 2D and 3D spectroscopy 
of H II regions has not revealed the existence of zones with sig- 
nificantly different abundances within the same H ii region, ex- 
cept in the case of NGC 5253, which shows some local N en- 
hancement (Walsh & Roy 1989; Kobulnicky et al. 1997, Lopez- 
Sanchez et al. 2007). On the other hand, chemical inhomo- 
geneities occuring on scales smaller than the spatial resolution 
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would pass unnoticed when comparing spectra of neighbouring 
zones. 

There is, however, a way to unravel the presence of small- 
scale abundance inhomogeneities in Hii regions. This is by 
comparing the abundances derived by traditional methods, i.e. 
from collisionally excited lines (CELs), to those obtained from 
optical recombination lines (ORLs) in the same spectrum. 
CELs are preferentially emitted in zones of high electron tem- 
perature and low metallicity, while ORLs are preferentially 
formed in zones of low electron temperature and high metallic- 
ity. So far, ORLs have been observed in only a few H ii regions. 
They systematically lead to an abundance discrepancy factor 
(ADFU larger than one. This has first been interpreted as due to 
temperature fluctuations within the H ii regions (Peimbert et al. 
1993; Peimbert 2003; Esteban et al. 1998, 1999a, 1999b, 2002, 
2004). However, Tsamis et al. (2003) consider that tempera- 
ture fluctuations alone cannot be the cause of the observed dis- 
crepancies, since in that case oxygen abundances derived from 



' The abundance disprepancy factor (ADF; Tsamis et al. 2004) of 
an ion is defined as the ratio of its abundance as derived from the 
intensities of ORLs to that derived from CELs. 
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far-infrared lines should be close to the high values obtained 
from optical recombination lines, and this is not the case (how- 
ever, aperture problems are difficult to deal with, see Garcia- 
Rojas et al. 2006). Tsamis et al. argue for the existence of a 
hitherto unseen component in H ii regions, consisting of cold, 
metal-rich ionized parcels of gas. In fact, chemical inhomo- 
geneities have been proposed for over a decade to explain the 
ORL/CEL discrepancy in planetary nebulae (Torres-Peimbert 
et al. 1990; Liu et al. 2000, 2004), but such an interpretation in 
the case of H ii regions has appeared only recently (Tsamis et 
al. 2003; Tsamis & Pequignot 2005). Of course, chemical in- 
homogeneities in H ii regions and planetary nebulae must have 
a very different origin. Note that invoking chemical inhomo- 
geneities accounts at the same time for the origin of tempera- 
ture fluctuations, which otherwise are difficult to explain quan- 
titatively (see Stasinska 2007 and references therein). 

In this article, we explore the possibility that the ADFs 
measured in H ii regions result, as suggested by Tsamis et al. 
(2003, 2005), from the presence of metal-rich "droplets'H as 
predicted in ffie scenario of Tenorio-Tagle (1996; hereafter re- 
ferred to as T-T96; see also Stasinska et al. 2007) for the en- 
richment of ffie ISM by Type II supernovae. Note that, up 
to now, reliable abundance discrepancy factors in H ii regions 
have been measured only for oxygen. For simplicity, we will 
ffien consider only oxygen in this paper Section 2 reviews the 
framework that accounts for the origin of such inhomogeneities 
in H n regions and evaluates their survival time-scale until full 
mixing with the ISM. This section gives also some observa- 
tional arguments in support of ffie model. Section 3 estimates 
how much matter can be expected in the form of metal-rich 
droplets. Section 4 discusses the ORL/CEL discrepancy in Hii 
regions in ffie context of metal-rich droplets. Section 5 presents 
some final comments and prospects. 

2. Metal-rich droplets in Hii regions 

2.1. Formation of the droplets 

The scenario proposed by T-T96 accounts for the fact that the 
ejecta from type II supernovae (SNe) ought to follow a long 
excursion into the galactic environment before they are able to 
mix wiffi ffie ISM. The excursion is promoted by the clustering 
of massive stars and thus of type II SNe, which, being corre- 
lated in space and time, force the violently ejected matter to 
generate large-scale superbubbles. These are able to exceed the 
thickness of galactic disks and burst into the haloes of the host 
galaxies, while displacing and locking the surrounding ISM 
into kpc-size expanding supershells. Such remnants, driven by 
ffie hot superbubble interior grow for as long as massive stars 
continue to release their metals, until the last 8 Mq star belong- 
ing to the star cluster or OB association completes its evolution 
(~ 40 Myr). During this time, the metals injected into super- 
bubbles have enough time to mix with the matter from stellar 
winds as well as with the matter thermally evaporated into the 

^ The word "droplets" is used in this paper with the same meaning 
as in Tenorio-Tagle (1996). It does not imply that they are in liquid 
form, but merely that they are of lower temperature and higher density 
than their surroundings. 



superbubbles from their suiTounding outer supershells. Mixing 
within the superbubble interior is strongly favored by the high 
temperatures (T ~ 10^-10^ K) and high sound speeds and by 
the stirring caused by the bursting of more SNe. However, the 
remaining ISM, being out of contact with the superbubble in- 
terior, is left uncontaminated by the products of the evolving 
star cluster After the last SN explosion, the low density within 
the superbubble interior (n ~ 10^^-10"^ cm"^) delays the im- 
pact of radiative cooling (fA ~ kT/{An) ~ 10^ yr; where A 
is the cooling rate and k ffie Boltzmann constant). Note also 
that radiative cooling does not occur at the same rate within 
the entire volume of the superbubble, but rather zones with a 
higher density will cool faster than their surroundings. Such a 
situation has been shown to lead, in other astronomical circum- 
stances such as galaxy formation or globular cluster formation 
(see Zel'dovich & Raizer 1966, Vietri & Pesci 1995), to mul- 
tiple re-pressurising shocks which in our scenario would lead 
to the formation of small, dense "cloudlets" of metal-rich gas, 
in pressure equilibrium wiffi the remaining hot gas. The high 
density of the cloudlets causes them to recombine, further re- 
ducing their temperature. The cold cloudlets would inevitably 
begin to fall towards the disk of the galaxy under the action 
of gravity. As ffiey fall through the lighter interstellar medium, 
they will be subject to various instabilities, such as Rayleigh- 
Taylor and Kelvin-Helmholtz (e.g., Rozyczka & Tenorio-Tagle 
1985; Schiano et al. 1995), which will lead to their breakup into 
a swarm of minute droplets. 

The overall outcome of such a fountain with a spray (as de- 
scribed by Scalo & Elmergreen 2004) is a large-scale dispersal 
of the matter processed by ffie star cluster over a large (kpc- 
scale) galactic volume, a volume much larger than the size of 
the stellar association. This is very different from other foun- 
tain models in which the hot gas rises to great altitudes (~ 
3 kpc) above the galactic plane to then cool all at once (see 
Kahn 1991). Such fountains would lead to massive giant clouds 
falling unimpeded towards the disk of galaxies. Note also that 
none of those models have addressed the mixing of metals with 
the ISM. On the other hand, the fountain with a spray of T-T96, 
predicts a stratus of dense, highly metallic and perhaps molec- 
ular droplets interspersed in the galactic ISM, wherever it had 
rained. Once this happens, the droplets are to participate in the 
general motion of clouds and intercloud medium promoted by 
galactic rotation and stellar activity, as envisaged by Roy & 
Kunth (1995). This is to lead to an even larger dispersal, al- 
though not to mixing of the metals with the ISM. For total mix- 
ing a new episode of stellar formation seems required as this 
will cause a thorough dissemination of the metals into the ISM 
enhancing its metallicity. This is primarily due to the UV radia- 
tion field, able to disrupt molecules and ionize the droplets and 
their surroundings, while generating multiple localized cham- 
pagne flows (Tenorio-Tagle 1979) and with them the stirring to 
promote a rapid mixing within the ionized volume. 

2.2. Survival of the droplets in H ii regions 

When the droplets are ionized, they will quickly (on order of 
a sound-crossing time) reestabhsh pressure balance with their 
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surroundings. The two processes that will tend to disperse the 
high concentration of metals in the droplets are molecular dif- 
fusion and turbulent mixing. 

The strong gradient in the concentration of oxygen ions be- 
tween the droplets and the surrounding gas will eventually be 
diffused away due to the collisional random walks of the parti- 
cles. The time required for diffusion through a distance L is of 
order 

fd = L^ID, (1) 

where D is the diffusion coefficient, which can be considered 
as the product of the root-mean-square particle velocity and the 
mean-free path between collisions (e.g., T-T96). For the diffu- 
sion of O^^ ions moving in a field of ions, D can be calcu- 
lated from the formulae presented in Oey (2003) as 

D = 1.04 X lO'^r^^^n"'/"' cm^ s"', (2) 

where T4 is the gas temperature in units of lO** K, n is the hy- 
drogen ion density, and / = 1 -i- 0.029 log Ti, — 0.010 log n. The 
resultant diffusion time is 

fd = 3.05 X lO^Lj^n/r;^^^ years, (3) 

where L\s is the length scale in units of 10'^ cm. Adopting 
characteristic values of density and temperature for O-rich 
droplets in pressure equilibrium (n = 1000 cm"^, T4 - 0.6, 
see Section 4), we find = lO^Ljj years. Therefore, metal-rich 
droplets as small as lO''* cm can survive over the entire lifetime 
of a typical H 11 region, but smaller droplets would be erased on 
smaller timescales and therefore would not contribute to the 
spectrum of any but the youngest H 11 regions. The correspond- 
ing mass of each droplet would be ~ 10 " M©. However, it 
is possible that disordered magnetic fields in the ionized gas 
might increase the diffusion time still further and allow smaller 
droplets to survive. 

Discounting this last proviso, we have a rather stringent 
condition on the properties of the droplets: a physical size of 
10^"* cm corresponds to an angular size of 0.015" at the dis- 
tance of the Orion nebula. Hence, droplets that are only a few 
times larger than this size may in principle be resolvable via 
direct imaging with the Hubble Space Telescope. There is al- 
ready some observational evidence that inhomogeneities may 
indeed be present at this scale (O'Dell et al. 2003) 

Diffusion by turbulent motions can be treated in a simi- 
lar way. Assuming turbulent eddies of velocity Hturb and size 
^turb, then for times less than the coherence or turnover time, 
^turb/fturb, matter is simply advected at a speed iimrb and so 
the effective diffusion coefficient grows linearly with time. For 
t > Auib/ftmb, on the other hand, the diffusion coefficient satu- 
rates at a constant value: 

i^turb - ^turbfturb- (4) 

However, this equation only gives the rate of dispersal of the 
metal-rich gas by the turbulent eddies (Klessen & Lin 2003), 
but does not directly address the physical mixing of the two 
phases (de Avillez & Mac Low 2002). Mere dispersal of the 
droplets is not enough to destroy their physical integrity and 
distinctive temperature. Instead, it is necessary for the turbulent 



eddies to stretch and distort the droplets, so as to steepen the 
concentration gradients until they reach a level at which molec- 
ular diffusion can occur (see above). The relevant timescale for 
this process can be derived from the exponential stretching of 
fluid elements by the turbulent shear (Pan & Scalo 2007) as 

tmix~td\nP. (5) 

In this equation, td is the turbulent diffusion timescale calcu- 
lated using equation ^ with ^turb = L (only those turbulent 
motions on the same scale as the droplets can efficiently dis- 
tort them) and !P - Dt^y\,/D is the Pec let number of the turbu- 
lence, with D given by equation (|2]l. If turbulence of velocity 
v' is driven at some large scale {', then the turbulent motions 
at a smaller scale ^tuA have velocity Oturb - ^'(Aurb/^')''. where 
fl = 1/3 for incompressible Kolmogorov turbulence (applica- 
ble to subsonic flows) and a - 1 /2 in the limit of highly super- 
sonic turbulence. The turbulence in H ii regions is likely driven 
by the mutual interactions of multiple transonic photoablation 
flows (Mellema et al. 2006; Henney 2006), with Mach numbers 
in the range 2-4 and at scales of 0.1 to 1 parsec. We therefore 
find a turbulent mixing time of 

W = 12, 600(Li5{[sy%l,g years, (6) 

where ^'jg and iiio are the scale and velocity of the turbulent 
driving in units of 10"^ cm and 10 km s"', respectively, and a - 
1/2 is assumedH This result implies that turbulent mixing is 
efficient on timescales shorter than those characteristic of H ii 
region evolution. 

Since the turbulent fd decreases as the droplet size is re- 
duced, the turbulence will quickly fragment the droplets to such 
an extent that the molecular diffusion also becomes efficient (at 
scales of order 10'^ cm), at which point chemical mixing will 
be complete. The long-term survival of the droplets in the Hn 
region therefore depends critically on the avoidance of turbu- 
lence. 

The extent to which turbulence is present in Hii regions 
is somewhat controversial. It is certainly the case that the in- 
ternal kinematics of H ii regions are highly irregular, with sig- 
nificantly non-thermal linewidths (O'Dell & Castaneda 1987). 
However, it is unclear whether this represents true hydrody- 
namical turbulence, or simply the superposition of multiple ve- 
locity components, each one of which arises in an ordered flow. 
Recent simulations of H ii region evolution in a clumpy molec- 
ular cloud (Mellema et al. 2006) indicate that the "chaotic" ve- 
locity field is largely the result of inward-pointing photoabla- 
tion flows, which form wherever the outward propagation of 
the ionization front is detained by a dense neutral condensa- 
tion. These flows shock against one another in the interior of the 
H II region, which does become truly turbulent. However, it is 
found that 10-30% of the [Oiii] emission comes from the non- 
turbulent flows closer to the ionization front. A similar result 
was derived analytically by Henney (2003), who showed that 
the surface brightness of such photoablation flows is propor- 
tional to the Unear size of the cavities that they carve out of the 

^ The expression also includes the term g = I + 0.034 log n^f + 
0.052 log Li5 - 0.017 log - 0.086 log 7+, but this will always be 
close to unity 
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H II region, which in turn is a function of the scale of the irregu- 
larities in the front. Flows from condensations with size > 10% 
the radius of the H ii region will have a brightness compara- 
ble to that of the entire region. Both these studies considered 
only H n regions with globally closed geometries, although it is 
likely that many optically visible regions are instead globally 
open champagne flows (e.g., Tenorio-Tagle 1979; Henney et al. 
2005), which would tend to increase the relative contribution of 
the photoevaporation flows to the total emission. On the other 
hand, the interaction of the stellar wind from the ionizing star 
with the other flows (e.g., Garcia-Arredondo et al. 2001) may 
be an additional source of turbulence, as may instabilities of 
the ionization front (Garcia-Segura & Franco 1996; Williams 
1999, 2002). Here we thus assume that turbulence and molecu- 
lar difiiision act very efifectively within the HII region volume, 
causing a thorough mixing and homogeinizing of the ionized 
gas. However, we also assume that at all times there are new 
metal-rich droplets traversing the ionization front into the HII 
region. 

2.3. Observational considerations 

There are several further observational facts that lend support to 
the droplet scenario, facts that ought to be taken into considera- 
tion by other possible explanations of how the mixing of heavy 
elements with the ISM takes place. Regarding abundances in 
Hn regions the issues are: 

Given the lifetime of type II SN progenitors (up to several 
10^ yr), it is obvious that the metals expected from recent bursts 
of steUar formation (age ~ a few 10^ yr) have not yet been 
released and thus the detected abundances in H n regions reflect 
the production, dispersal and mixing of heavy elements ejected 
by former stellar generations. 

This result is supported by the work of Herrero & Najarro 
(2005) and Simon-Diaz et al. (2006) who measured the metal- 
Ucity of the ionizing stars (in Orion and in the largest H n re- 
gions in M33) and found that these are similar to those deter- 
mined for their corresponding H n regions. Thus the material 
gathered by the star(s) during collapse and the star formation 
process had the same contaminating metals and in the same 
proportion - given by Eq. (22) - as those detected in their pho- 
toionized volumes. 

Low mass galaxies such as I Zw 18 (Vflchez & Iglesias- 
Paramo, 1998; Legrand, 2000), II Zw 40, NGC 4214 
(Kobulnicky & Skillman, 1996), NGC 1569 (Devost et al. 
1997), Sextans A and B (Kniazev et al. 2005), in which abun- 
dances have been measured in several locations all present very 
uniform metal abundances within their ionized volumes. Given 
the size of these H n galaxies (> 1 kpc), these facts imply a very 
even large-scale dispersal of the metals from former stellar gen- 
erations. Scalo & Elmegreen (2004) summarize the evidence 
for abundance variations in the ISM of our Galaxy, concluding 
that the gas is well-mixed on scales from 1 to 100 parsec. Note 
that the inhomogeneities posited in the scenario considered in 
this paper are on much smaller scales. 

Abundances in the ISM can also be obtained from the ab- 
sorption lines produced by the Hi envelopes that surround Hn 



regions. These derivations are difficult and uncertain and re- 
quire very high resolution far-UV spectroscopy. The more reU- 
able, recent results in external galaxies (e.g. Aloisi et al., 2003, 
Lebouteillier et al. 2004) indicate an oxygen abundance in the 
H I regions several times lower than in the H n regions (although 
this is not always the case, see Thuan et al. 2005). A systematic 
difference between abundances in H n and H i regions has many 
possible explanations. One of them is that abundance determi- 
nations using absorption and emission lines have significantly 
different biases in the case of abundance inhomogeneities. 

Finally, using long-slit spectroscopy of the Orion nebula, 
Rubin et al. (2003) have shown the presence of temperature 
fluctuations on a scale of 0.5". In the same object, from the ra- 
tio of [Om] /14363 and [Om] ^007 Hubble Space Telescope 
images, O'Dell et al. (2003) found small scale temperature fluc- 
tuations of + 400 K on the 10" scale, but with most power be- 
ing at smallest scales: 0.1", i.e., 6.45 x lO^^^cm. O'DeU et al. 
(2003) say that shadows behind clumps close to the ionization 
fronts can partly account for the observed temperature fluctua- 
tion, but an additional cause is needed. As a matter of fact, these 
observed fluctuations could be due to our metal-rich droplets. 
The dominant scale of the observed fluctuations is tantalisingly 
close to the droplet sizes derived in Section 2.2. Correlations of 
these observed fluctuations with surface brightness would give 
a hint on whether our interpretation is correct. 

3. How much does it rain? 

In order to see what mass fraction of oxygen can be expected 
in the droplets, here we consider a simple scheme for the life 
cycle of metal-rich droplets assuming a steady state. 

3.1. Formal derivation 

If one assumes that all the oxygen ejected by massive stars 
eventually rains down in the form of highly metallic droplets, 
the rainfall rate of oxygen in drops at a given time t is given by: 

R(t) = e.(0) i,(t - At), (7) 

where e,(0) is the mass fraction of all stars formed that is 
ejected as oxygen by SNe, and S»(f - At) is the star forma- 
tion rate (SFR) per unit area at t - At, with Af ~ 10^ yr, the 
time needed for SNe ejecta to travel through the halo and re- 
turn to the disk. The droplets end up being homogeneously dis- 
tributed in the ISM, to then be fully processed while enhancing 
the metalhcity of newly formed stars and of their photoionized 
surrounding gas, every time that they participate in new bursts 
of steUar formation. 

If Xi is the mass of gas that is ionized per unit mass of stars 
formed, the fraction of gas ionized or incorporated into stars 
per unit time is lit) = (Xi + 1) S,(?)/2ism(0. where Sjsm is the 
surface density of gas. The rate per unit area at which oxygen 
is released from the droplets to the ISM is: 

m 2d(0; t) = ^ Ed(0; t), (8) 
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where £d(0; t) is the oxygen mass in droplets per unit area at 
time t. The time variation of 2d(0; t) is then: 



di:D(0; 



(9) 



If one considers a time interval t - to during which 2t(f - At), 
2,(0, and SismCO are approximately constant, eq. (9) can be 
integrated to give: 



lD(0;0 = (sD(0;fo)-y) 



The drops processing timescale is thus: 

_ 1_ _ ^ISM 



(10) 



(11) 



If we take Xi - 9 masses of ionized gas per unit mass in new 
stars (Williams & McKee 1997), with Sjsm - 10 Mq pc"^ 
(see Fig. 16 of Molla & Diaz 2005) and t, ^ 1.5 Mq pc'^ 
Gyr ' (McKee & Williams 1997) in the solar neighbourhood, 
td - 130 Myr. And for t-to equal to 2-3 times the processing 
timescale td, we get from eq. (10): 



V A ^ e.(O)2.(f-A0SisM 
ED(0;n ^ — = : , 

/ (Xi + l)S,(0 



(12) 



and if the SFR has been approximately constant for the last one 
or few Gyr (probably a good approximation to within a factor 
of 2 for our own Galaxy - see Fig. 15 of Molla & Diaz 2005): 



^D(0;t) 



e,(0) £isM 
(^, + 1) ■ 



(13) 



The surface density of gas in the droplets with respect to 
the surface density of gas in the ISM is 



X = 



2d 



SD(0;f) Zism(O) 



SisM Sism(O;0 Zd(0) 



(14) 



where Zd(0) and Zism(O) are the oxygen mass fractions in the 
droplets and in the ISM, respectively. 
If we define the parameter rj by: 



n = 



(Xi + l)Z0(O)' 
it follows from Eqs. (13) and (14) that: 

n 



X = 



Zd(0)/Zo(0)- 



(15) 



(16) 



In first approximation, the value of rj is not expected to vary 
strongly with metalUcity, so the value of X is roughly inversely 
proportional to the metallicity of the droplets. 

3.2. Numerical estimates 

Neither the value of e,(0) nor the value of Xi are well-known. 
In Table 1, we show different estimations of e,(0). This is ap- 
proximately equal to the mass fraction of oxygen ejected by 
stars of all masses, since the contribution of intermediate mass 
stars is negUgible. The values shown in Table 1 are based on 
different sets of calculations of stellar yields and on different 



Table 1. Mass fraction of all stars formed that is ejected as 
oxygen. 



TA /Tea 
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Z/Zq 


e,(V) 


Kei 


Scalo 
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Scalo 


1 
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O.UOzz-U.OOz / 
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Scalo 


U.Uz 


U.U04U 
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Scalo 
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A AA1 A 
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2 


Salpeter 


n 1 
U.i 


A AA^T A AAA1 

U.UUjZ— U.UUuZ 


1 


Salpeter 


1 


0.0056-0.0070 


1 


Salpeter 


0.02 


0.0104 


2 


Salpeter 


1 


0.0100 


2 


Kroupa 


0.1 


0.0092-0.0110 


1 


Kroupa 


1 


0.0099-0.0120 


1 


Kroupa 


0.02 


0.0184 


2 


Kroupa 


1 


0.0177 


2 



1- Woosley & Weaver (1995); 2- Portinari et al. (1998) 
^ The IMF slopes are those of Scalo (1986) - as given in 
Lanfranchi & Matteucci (2003) -, Salpeter (1955), and 
Kroupa (2001) 



stellar initial mass functions. There is presently no consensus 
on which are the most realistic values and chemical evolution 
models of galaxies have been constructed exploring different 
sets of yields and initial mass functions (e.g., Portinari et al. 
1998, Henry et al. 2000, Chiappini et al. 2003, Francois et al. 
2004). In Table 1, the values of e,(0) range between 0.002 and 
0.018. They have been computed for a lower stellar mass of 
Mdown - OA Mq and an upper stellar mass of M^p - 120 Mq. 
Adopting larger values of Mdown but maintaining the same ini- 
tial mass functions would lead to somewhat larger values of 
e,(0) (e.g. 0.015 for Mjown = 0.3 Mq instead of 0.01 for 
Mdom =0.1 Mq in the case of the yields of Portinari et ill. 
1998 with the Salpeter initial mass function). 

Concerning Xi, the situation is probably even more com- 
plex, because an estimate of this parameter involves a repre- 
sentation of the structure of the ISM on different scales and 
a knowledge of all the forces at play, in addition to the IMF 
description. Using a simple argument of ionization balance, 
Izotov et al. (2006) derive Xi-360/n, where n is the ionized gas 
density. So, for a density of 300 cm"^, Xi would be close to 1. 
In a comprehensive study of the evolution of molecular clouds 
in the Galaxy, Wilhams & McKee (1997) give Xi = 9 while 
from the work of Franco et al. (1994) one derives Xi ~ 25. The 
parameter Xi is likely to change with metallicity due to a vary- 
ing amount of ionizing photons absorbed by dust or to a change 
in the hardness of the ionizing radiation field. 

Figure 1 shows the variations of X as a function of Zd for 
T] = 0.17, our working value in the remaining of the paper 
(which corresponds to e,(0) = 0.01 and Xi = 9). 

Finally, note that a change by a given factor in the star for- 
mation rate between t - At and t would change X by the same 
factor. 
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Table 2. Densities, temperatures, metallicities and ADFs measured in H ii regions 



ID 




re[Oiii](K) 


TeCBJ) (K) 


12 + logCO/H)" 




ADF(0++) 


Ref 






Galactic H ii leg 


ions 










M8 


1750 


8050 ± 700 




8.41 


0.52 


2.0 


1 




1800 


8090 ± 140 


7100""'^'° 

' ^""-1000 


8.51 


0.66 


2.3 


2 


M16 


1120 


7650 ± 250 


5450 ± 820 


8.56 


0.74 


2.8 


3 


M17 


860,520 8120 ±250, 8210 ±250 




8.53, 8.50 


0.69, 0.65 


1.8, 2.2 


4 




600-1500 


8200 


7700 


8.56 


0.74 


2.1 


5 




470 


8020 ± 170 




8.52 


0.68 


1.9 


2 


M20 


270 


7800 ± 300 


6000 ± 1300 


8.53 


0.69 


2.1 


3 


M42 


4000, 5700 8300 ± 210, 8350 ± 200 8730 ± 800, 8390 ± 800 


8.47, 8.47 


0.60, 0.60 


1.3, 1.5 


6 










8.52 


0.68 


1.3 


5 




8900 


8300 ± 40 


7900 ± 600 


8.51 


0.66 


1.4 


7 


NGC 3576 


1300-2700 


8850 


8070 


8.52 


0.68 


1.8 


5 




2800 


8500 ± 50 


6650 ± 750 


8.56 


0.74 


1.9 


8 


NGC 3603 


5150 


9060 ± 200 




8.46 


0.59 


1.9 


3 


S311 


310 


9000 ± 200 


9500 ± 900 


8.39 


0.50 


1.9 


9 


Extragalactic H ii regions 


LMC 30 Dor 


370-1800 


10100 




8.34 


0.45 


2.0-2.7" 


5 




316 


9950 ± 60 


9220 ± 350 


8.33 


0.44 


1.6 


10 


LMCNllB 


80-1700 


9400 




8.41 


0.52 


4.9-8.2'' 


5 


NGC 604 


< 100 


8150 ± 150 




8.49 


0.63 


1.6 


11 


NGC 2363 


360-1200 


15700 ± 300 




7.87 


0.15 


2.2 


11 


NGC 5253 


370-610 


10940-12010 




8.18-8.28 


0.31-0.39 


1.5-1.9 


12 


NGC 5461 


300 


8600 ± 250 




8.56 


0.74 


1.9 


11 


NGC 5471 


220-1150 


14100 ± 300 




8.03 


0.22 


1.6 


11 


NGC 6822 V 


175 


11900 ±250 




8.08 


0.24 


1.9 


13 


SMC N66 


50-3700 


12400 




8.11 


0.26 


2.3 


5 



1- Esteban et al. (1999a); 2- Garci'a-Rojas et al. (2006b); 3- Garci'a-Rojas et al. (2006a); 4- Esteban et al. (1999b); 

5- Tsamis et al. (2003); 6- Esteban et al. (1998); 7- Esteban et al. (2004); 8- Garci'a-Rojas et al. (2004); 

9- Garcfa-Rojas et al. (2005); 10- Peimbert (2003); 11- Esteban et al. (2002); 12- Lopez-Sanchez et al. (2007); 



13-Peimbertetal. (2005) 

" The high values of the densities in the extragalactic Hn regions mostly come from the [Ariv] ratio and are very uncertain 

12 + log(0/H) from forbidden lines 
° the Solar oxygen abundance is taken from Allende Prieto et al. (2001) 

■i In 30 Dor and Nl IB, the highest values of the ADFs include a correction for contamination by scattered light 




Fig. 1. The value of X = Md/Mism as a function of the metal- 
licity of the droplets, assuming that tj = 0.17 (see Eqs. 18 and 
19 in text). 



4. Metal-rich droplets and the ORL/CEL 
discrepancy 

4.1. ADFs measured in H ii regions 

Table 2 lists the values of the densities, temperatures and abun- 
dance discrepancy factors ADF(0^^) found in the literature 
for Galactic and extragalactic Hii regions. The temperatures 
derived from Balmer jump, re(BJ), are not available for all 
cases and bear generally significant uncertainties. It is however 
clear that they are, on average, smaller than those derived from 
[Om] A4363/5007 line ratios, T^[Om]. The hsted ADFs are 
equal to the ratio of the value of the O"^"^ abundance derived 
from [O III] /15007 to the one derived from O ii recombination 
lines, using T^[0 III] in both cases. These ADFs were obtained 
from several recombination lines of O ii and the forbidden line 
[Oiii] /15007, using the temperature rr[0 iii] derived from the 
[Om] ^4363/5007 ratio. The values of ADFCO^-') are found 
to be around 2, over a large range of H n region metallicities: 
12+log(0/H) = 7.87-8.56 (i.e., Z = 0.15 to O.TSZq). However, 
note that the ADFs measured for most of the extragalactic Hn 
regions in Table 1 and, in particular, for all the low metallicity 
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objects (12 + log (O/H) < 8.3) were derived using the O ii lines 
of multiplet 1, and these Unes can be severely affected by stel- 
lar absorption lines (Tsamis et al. 2003). This means that the 
ADFs could be much larger than 2 for these objects. 

4.2. A simple expression for tlie ADF 

To better understand what an ADF really is, we derive here an 
approximate expression for it in the case of a two-component 
toy model. In a medium characterized by an electron density We 
and an electron temperature T^, the luminosity in an optically 
thin line arising from O"^"^ ions is given by 



(17) 



where NiO'*'^) is the total number of O"^"^ ions in the emitting 
volume, and e{Te) is the line emissivity, in erg cm^ s"\ which 
depends essentially on the electron temperature (for the sake 
of simplicity, we neglect the density dependence, as it is not 
important in the problem considered here). The expression for 
ADF(0^^) is then given by 

^ORl/ fORL(7'e) 



ADF(0'^"') = 



(18) 



where is taken equal to r, |0 iii]. However, if within the 
emitting volume there are two different media with different 
physical parameters such as their metalUcity, then 



ADF(0"'"') 



ORL 



'CEL 



+ -L^EL)/W7'e)' 



(19) 



where a and b refer to the two media, and is here again equal 
to the temperature derived from the observations, Tj{0 iii]. In 
the limiting case in which the medium b is too cold to con- 
tribute significantly to the CEL intensities (which is easy to 
achieve due to the very strong temperature dependence of the 
/i4363 and A5007 lines ), Tr\0 m] is in fact the temperature of 
medium a. Using Eq. (17), Eq. (19) becomes: 



ADF(0"'+) = 1 + 



or. 



ADF(0"'"') = 1 -I- 



A^\0^n|eoRL(7t) 
iV«(0++)<eoRL(re")' 



AO^n|Z*M^eoRL(7t) 
x«(0++) nlZ" M- eoRi.{T^y 



(20) 



(21) 



where x", yf' stand for the doubly ionized oxygen fractions, 
Z", are the O/H ratios (or "metallicites"), and M", are 
the total hydrogen masses in zones a and b, respectively. Note 
that eoRL is an increasing function of Z via T^, and x{0^^) is a 
decreasing function of ne due to recombination. The observed 
values of the ADFs in H ii regions impose strong constraints on 
the parameters that come into play. 

4.3. Pliotoionization modelling 

We now simulate the situation described in Sect. 2 by a set 
of photoionization models. We consider three ionized com- 
ponents, each with its own metalhcity: the highly metalUc 
droplets, the background ISM, which is being ionized at the 




Fig. 2. The metallicity of the fully mixed H ii region as a func- 
tion of the metallicity of the droplets, for Zism/Z© = 0.2 (cir- 
cles) and Zism/Z© = 0.6 (squares), assuming that t; = 0.17 (see 
Sect. 3.1). 

same time as the droplets, and the H ii region into which the 
former two media have mixed. At each moment, a mass Mq 
of droplets with metallicity Zd and a mass Mjsm from the ISM 
with metallicity Zism are being ionized and will eventually re- 
sult in a mass Mmi of ionized gas with metalhcity 

XZd + ZjsM 



X-l- 1 



(22) 



where X = Md/Mism- Figure 2 shows Znn as a function of Zd 
for Zism = 0.2 and 0.6 Zq, taking 77 = 0. 17. 

The real progression of the ionization front is complex, but 
the situation can be modelled in a simple way by consider- 
ing that the ISM, the droplets and the fully mixed H 11 region, 
with typical densities hism^ and nnib respectively, are im- 
pacted by the same ionizing radiation field. This implies that 
the droplets do not have too large a covering factor so that the 
general ISM can be ionized (in principle, a situation where the 
droplets would completely block the ionizing photons is possi- 
ble as well, but not considered in this paper). 

Here we assume that any observation of an H 11 region en- 
compasses zones close to the ionization front that have just 
been ionized and still have their initial metallicities Zd and 
Zism, and zones that have been ionized some time ago and have 
acquired, through the turbulent stirring promoted by well local- 
ized champagne flows (see Sect. 2), the metallicity Zhu given 
by Eq. (22). We assume that this fully mixed gas contributes a 
fraction, F, of the total observed H/3 flux from the H 11 region. 
Note that the parameter F also accounts for the lifetime of the 
droplets once they are ionized. 

For a given stellar radiation field, and a given value of 77 - 
as defined by Eq. (15) - the problem is entirely specified by 
the parameters F, Zd and Zism, and the densities of the three 
different media. We expect the densities nisM and nun to be 
of the same order of magnitude because the fully mixed ion- 
ized gas must be in rough pressure equihbrium with the ISM 
gas that has just been ionized and because temperature differ- 
ences between the media are likely smaU, which is confirmed 
by our results shown in Figs. 4d and 5d. For the computations, 
inspired by the observed values listed in Table 2, we have taken 
«iSM = «Hn = 300 cm"^. Since the droplets are more metaUic 
than the ISM, their temperature must be lower, and thus we 
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Fig. 3. The average temperature in the O^^ zone of the droplets, 
as a function of the metalhcity of the droplets. 



assumed the droplets to be denser than the background. We 
have considered various density contrasts and present results 
from calculations assuming - 1000 cm"^. The computa- 
tions were done with the code PHOTO (Stasiriska 2005), taking 
for the ionizing radiation field that of a blackbody at a tempera- 
ture TefF = 50,000 K. The inner radius of the models is the same 
for the three components, and has been chosen so that the ion- 
ization parameter at the inner radius, U, is equal to 10"^ both 
for the H ii region and the Fig. 3 shows how the O^^- 

weighted temperature in the droplets decreases as a function of 
the droplet metallicity. 

We computed models of the droplets with various values of 
Zd, and models for the mixed H ii region with Zhh as given by 
Eq. (22), the value of X resulting from Eq. (16) for an adopted 
J] of 0.17. We combined the photoionization models for the 
droplets, the ISM and the fully mixed Hii region taking var- 
ious values of F. We then applied to the resultant models the 
same techniques as an observer deriving Tj[0 m], ADF(0^^), 
and the metallicities from the coUisionally excited lines and re- 
combination lines. We used exactly the same atomic data as for 
the photoionization models. 

4.4. Discussion of the simulated observations 

The number of parameters in this problem is large, and dis- 
playing our results for the entire parameter space would be 
cumbersome. Since the main goal of this paper is to investigate 
whether the ADFs observed in H ii regions could indeed be the 
signature of the enrichment scenario proposed by T-T96, we 
will concentrate on a few plausible cases, and try to learn some 
lessons from them. 

Figure 4 presents the results of simulations for Zism = 
0.2 Zq, taking 77=0.17 (the working value presented in Sect. 
3.2). In all the nine panels the abscissa represents Zq in solar 
units. The range of values for Zd was inspired by the work of 
Silich et al. (2001) where the metallicity of the hot superbubble 
interior, resultant from the multiple SN explosions of an aging 
star cluster, was calculated. The various cases they considered 
(see their Fig. 6) encompass starbursts evolving in low and high 



metallicity galaxies and account for mass evaporation from the 
surrounding supershell into the hot superbubble interior. The fi- 
nal values, the metallicity of the matter that will eventually cool 
to compose the high metallicity droplets, range from 0.2 Zq to 
5 Zq. In the plots, simulations corresponding to an observation 
where F = 0.1 are represented by circles, those corresponding 
to F - 0.5 by squares, and those corresponding to F - 0.9 
by triangles. Panel a shows the simulated ADF(0^^), panel b 
shows Ti[0 III], panel c shows Zqel, the metallicity derived 
from coUisionally excited lines. These three quantities repre- 
sent "observables". The next panels show other results from the 
simulations, which shed some light on the behaviour of the "ob- 
servables". Panel d shows r(0^^)Hii, the O^^-weighted tem- 
perature in the fully mixed H n region (which, as can be seen 
by comparing panels b and d, is not identical to the temper- 
ature derived from the [Om] /14363/5007 ratio). Panels e and 
/ show the bias in the determination of the metallicity when 
using coUisionally excited lines (Zcel/Zhii) or recombination 
lines (Zorl/Zhii). The computed values of Zcel and Zqrl take 
into account the contribution of O^, which has been evaluated 
from the [O ii] A3121 line using the temperature derived from 
the [Nil] /15755/6584 line ratio. Note that this procedure gives 
lower ORL abundances than the procedure used by most au- 
thors quoted in Table 2, who correct the temperature of the O^ 
zone for temperature fluctuations using the scheme of Peimbert 
(1967). Panels g, h and / show the flux from the droplets in the 
[Oiii] /i4363, [Oiii] A5QQ1 and On /1465 1 El Hnes, with respect 
to the total in these lines. 

As expected, ADF(0^^) increases with Zd, and, for a given 
Zd, it is larger for smaller values of F. The computed Zcel also 
strongly depend on F. In other words, the computed metal- 
licity is higher if the observation gathers less light from the 
fully mixed H ii region with respect to the newly ionized gas. 
The behaviour of 7'r[0 iii] is more intricate. As Zd increases, 
the temperature in the drops becomes smaller, as seen in Fig. 
3, due to the increased cooling by metals. One would expect 
7'r[0 III] to decrease as Zd increases. This is seen only at small 
values of Zd. The behaviour of rr[0 iii] is actually due to a 
subtle combination of various effects, linked to the tempera- 
tures in the three phases of our models and in the respective 
proportions of these phases in the resultant model. The temper- 
ature in the fully mixed H ii region decreases slightly as Zd in- 
creases, because Zhu increases as prescribed by Eq. (22), which 
accounts for mass and metal conservation. However Fig. Ag 
shows that the drops contribute less to the [Oiii] /t4363 flux 
than to the [Oiii] /15007 flux, and thus, the [Om] ^4363/5007 
ratio is larger when the contribution of the resultant H ii region 
to [Om] /15007 is smaller. The value of Zcel is not only re- 
lated to the ADF, but also to the contribution of the O^ zone. 
Therefore, the behaviour of Zcel with Zd seen in Fig. 4 c is not 
necessarily universal. A wide range of numerical models, with 
a better representation of the physical model and a better sam- 
pling of excitations conditions would be needed to draw some 
useful conclusions from the computed values of Zcel- 



We use the definition U = Qn/(4nR^n), where is the number 
of stellar ionizing photons, R is the distance to the star, and n is the 
gas density. 



^ Following the usual nomenclature for recombination lines, the O n 
/14651 line is emitted by the recombining O^^ ion. 
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Fig. 4. Results of simulations for a series of models with Zism - 0.2 Zq as a function of Zd, for various values of F and taking 
?7=0.17 (see Sect.3.2 and 4.3). Circles: F=0.1; squares: F = 0.5; triangles: F = 0.9. "Observables": (a) ADF(0++), (b) T,[0 iii], 
(c) ZcEL (see text). Other results of the simulations: (d) r(0^^)Hii, (e) Zcel/Zhii, i.e., the metallicity bias when using [Oiii] /15007 
(see text) (f) Zorl/Zhii, i.e., the metallicity bias when using the O u /1 465 1 line to measure the oxygen abundance (see text); (g), 
(h) and (i): the flux from the droplets in the [Om] /14363, [Om] /15007 and On /14651 lines, with respect to the total in these 
lines. 



Figure 5 is the same as Fig. 4, but now for simulations with 
ZisM=0.6 Zq. Qualitatively, it displays the same characteristics 
as Fig. 4, but with less extreme features. The ADFs(0^'^) are 
here much smaller, because the temperature contrast between 
the metal-rich drops and the ISM is smaller. On the other hand, 
simulations with Zism smaller than 0.2 Z©, not shown here, lead 
to more extreme characteristics than displayed in Fig. 4. 

The same set of figures for droplets with a density riu = 
3000 instead of 1000 cm"^ are very similar to Figs. 4 and 5, but 
with smaller effects, due to the fact that gas is more recombined 



in the droplets than in the cases shown here. However, the den- 
sity effect is only slight, since, as seen in Eq. (21) the effects of 
density and recombination on ADF(0"^'^) tend to compensate. 

Our computations were carried out with a given stellar ra- 
diation field and with a given ionization parameter. We do not 
expect large changes in our results when modifying the steUar 
radiation field or ionization parameter within the limits rele- 
vant to bright Hn regions. The average effective temperature 
does not change much from one object to another, and is a sec- 
ond order parameter in our problem. The ionization parameter 
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Fig. 5. Same as Fig. 2 but for models with Zism = 0.6 Zq. 
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can be different from object to object, but the ADFs are abun- 
dance ratios of the same ion, and thus are not dependent on 
the ionization parameter to first order. Test computations have 
confirmed this view. 



The most important parameters are Zp and F (which are 
varied in the plots presented in Figs. 4 and 5) and our parameter 
7], which, from considerations in Sect. 3.2, has a large range of 
possible values. 



Tkble 3. Expected range of parameters entering the model. 



parameter 


values 


argumentation 


e,{0) 


0.002-0.018 


Table 1 


Xi 


1 ? - 25 ? 


see Sect. 3.2 


Zd 


Zism - lOZ©? 


Silichetal. (2001) 




100 -5 000 cm-' 


Table 2 


"ISM /"HIl 


^ 1 


see Sect. 4.3 


"D MsM 


3-30 


rough pressure equilibrium 


u 


10-3 - 10-1 


Stasinska & Leitherer (1996) 


Tes 


35,000 - 55,000 K 


massive stars 


F 


0-1 


see Sect. 4.4 



Table 3 summarizes our knowledge on the range of param- 
eters entering the model. 
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4.5. What the confrontation between computed ADFs 
and observed ADF tells us 

Let us now compare the simulations shown in Figs. 4 and 5 with 
the observations. The Galactic and Magellanic Clouds nebulae 
Usted in Table 2 have diameters of (10-50)' while the shts have 
widths of (1 - 3)", so that the spectroscopic observations cover 
a tiny fraction of the objects, and the contribution from the fully 
mixed H n region is hkely small. The remaining extragalactic 
objects are smaller, with angular radii of (10 - 100)", so the 
contribution of the fully mixed H n region is more important, 
however still small. We thus believe that the models to be com- 
pared to the available observations are those corresponding to 
the smallest values of F. The ADFs listed in Table 2 are mostly 
in the range 1.5-2, with only few cases above 2. Such values 
are obtained for our models for a variety of physical conditions. 

Our models show clear trends in the values of ADFs with 
ZisM, while there is no obvious trend of the observed ADFs 
with any observed physical parameter. Note, however, that, for 
the extragalactic H n regions, the uncertainties in the absorp- 
tion corrections (and whether they should be applied) and the 
fact that for most of the extragalactic objects the ADF has been 
determined from just a couple of hues (or line-blends) of mul- 
tiplet 1 (which are very weak and difficult to measure), pre- 
clude us from excluding any of the following possibilities: i) 
all the extragalactic H n regions show ADFs ~ 2; ii) some ex- 
tragalactic H n regions have high ADFs; Hi) the low-metalhcity 
H n regions show higher (lower) ADFs. From our models, low 
metalUcity objects are bound to show higher ADFs. Obviously, 
a more conclusive sample of observational data is needed to be 
confronted with our model predictions. 

It has been commented by Lopez-Sanchez et al. (2007) that 
observed ADFs are similar in a variety of H n regions in differ- 
ent host galaxies. They take this fact as an argument against 
metal-rich droplets producing the ADFs. The entire grid of 
our models does predict a wide range of ADF values, which 
seems in contradiction with the very restricted range of ob- 
served ADFs and the fact that there seems to be no correla- 
tion between them and any observed property of the H n re- 
gions (see Garcia-Rojas & Esteban 2006). However, not all our 
models necessarily represent a physically plausible model. For 
example, the value of Zd is certainly constrained by physical 
processes occuring in the superbubbles, so that only a fraction 
of the considered range is relevant. In this respect, we might 
even suggest that, in the framework of our scenario to explain 
the ADFs, the observed values could put constraints on the 
mixing of metals in the supernova bubbles. For example, the 
observed ADFs do not seem to depend on metalhcity whereas 
our models would predict larger ADFs for a lower Zism if all 
other parameters were to remain fixed. From Table 2, the values 
of e,(0) do not change with metalhcity. Thus, if our scenario 
is to explain the observed ADFs, this imples that one or more 
of the remaining parameters (F, X;, or Zd) must correlate with 
Zism in such a way as to maintain a constant ADF. Indeed, one 
might expect a relation between Xi and Zism. via metallicity- 
dependent stellar winds which play an important role in the 
shaping of the ISM and therefore on determining the mass of 
gas that is ionized. One might also expect a relation between Zd 



and Zism, although Silich et al. (2001) show that this relation is 
not a simple one. 

Note that the magnitude of the biases found in our simu- 
lations obviously depend on the value adopted for rj. Had we 
used a value of ?7=0.013, a minimum value impUed by the ref- 
erences listed in Sect. 3.2, the biases would have been very 
small in most cases. 

4.6. Possible effects of dust 

In our models we did not consider the effect of dust. This 
was on purpose, in order not to inflate the number of free 
parameters. One of the expected effects of dust would be to 
decrease the amount of O available to produce recombination 
Unes in droplets if a significant fraction of the O atoms in drops 
are depleted onto refractory dust grains. However, maximum 
dust condensation efficiencies in SNe ejecta seem to be ~ 0.2 
(Todini & Ferrara 2001; Edmunds 2001), so the effective value 
of X(0), would only decrease to by a factor of 0.8. We be- 
heve that the main effects of dust grains in this problem is their 
contribution to heat the gas. The efficiency of grain heating de- 
pends essentially on the dust-to-gas mass ratio, the ionization 
parameter, and the size of the grains. The contribution of grains 
to the total heating budget is larger for higher ionization param- 
eters (Stasinska & Szczerba 2001). Since the droplets are to be 
denser than the surrounding gas, the ionization parameter there 
is lower, and therefore the effect of grain heating is expected 
to be smaller than in the surrounding gas. Thus, the effect of 
the ionization parameter alone would be to enhance the con- 
trasts between droplets and the rest, and thus to enhance the 
effects discussed in the preceding section. However, the prop- 
erties of the dust grains are likely very different in the droplets 
and in the interstellar medium, since the dust in the droplets 
should be heavily weighted by the conditions during the super- 
nova ejection and droplet condensation. This issue obviously 
merits further study. 

4.7. The abundance bias 

It is often said that abundances derived from recombination 
lines are not biased, because all recombination lines have 
roughly the same temperature-dependence (e.g., Esteban et al. 
2005). On the other hand, abundances derived from collision- 
ally excited fines could be underestimated if large temperature 
fluctuations in the O"^"^ zone are present but not accounted for. 
Some authors apply various techniques to try to estimate tem- 
perature fluctuations in a given object and correct for them us- 
ing the scheme of Peimbert (1967). Others do not apply such 
a correction, considering that temperature fluctuations are ex- 
pected to be small. 

In the case of our scenario, the real metallicity is that of 
the fully mixed H n region. Figures 4 e and 5 e show that the 
bias Zcel/Zhh is larger than 1 in all the cases shown here. 
Obviously, in our scenario, one cannot expect the abundance 
measured from collisionaUy excited lines to measure the real 
metalUcity exactly. The discussion of the abundance bias is un- 
fortunately complicated by the fact that one should also con- 
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sider the zone, and that the results are expected to strongly 
depend on the excitation of the object. We note that, in our 
models, the bias Zcel/Zhii is usually small, but it may reach 
0.2 - 0.3 dex if the weight of the well-mixed H ii region is small 
and if is of the order of 2 Z©. It is interesting to note that this 
bias is opposite to the bias due to the temperature fluctuations 
postulated by Peimbert (1967). However, many more numeri- 
cal experiments would be needed to consider this as a general 
rule. 

Figures 4/ and 5/ show that the bias Zorl/Zhii is always 
larger than one, and may reach quite large factors. As expected, 
the bias increases with ADF(0^^). 

To summarize, if the droplet scenario is valid, ORLs over- 
estimate the real metallicity of the ISM, possibly by quite 
large amounts. CELs may also overestimate it, although to a 
much weaker extent. The biases are stronger for larger ADFs. 
Unfortunately, we have at present no recipe to correct for the 
biases. In the meantime, a prudent attitude would be to discard 
objects with large ADFs as probes for chemical evolution stud- 
ies of galaxies. 

5. Conclusions and future directions 

In this paper, we have described a scenario for the enrichment 
of the ISM and discussed its observational signatures. The sce- 
nario invokes metal-rich droplets produced by a long process 
following supernova explosions in previous stellar generations. 
The full enrichment of the ISM is only achieved when these 
droplets are destroyed by diffusion in photoionized regions. 
Following Larson (1996), if the total star formation rate in the 
galaxy is about 3 Mq yr"' and the amount of cloud matter 
that can be ionized per every solar mass converted into stars 
amounts to ~ 10 Mq , then if the total gaseous mass of the 
galaxy amounts to 4 x 10^ Mq, the ionization rate implies that 
the entire ISM is cycled through the ionized phase every 100 
Myr. It is then within such a short timescale that the ISM ac- 
quires, at least locally, an even metallicity. On the other hand, 
the time required to disperse the metals from a stellar genera- 
tion is much longer (several xlO^ yr, due mostly to a slow cool- 
ing process) and a significant fraction of the mass converted 
into stars in the last few stellar generations, which ends up be- 
ing violently reinserted into the ISM, remains hidden within 
the hot phase of the ISM. 

We have shown in Section|4]that photoionization of highly 
metallic droplets can, under certain conditions, reproduce the 
observed abundance discrepancy factors (ADFs) derived for 
Galactic and extragalactic Hii regions. In this scenario the re- 
combination lines arise from the highly metallic droplets that 
have not yet fully mixed with the ISM gas. However, we find in 
Section lZSl that the droplets in H n regions can only have a nar- 
row range of sizes, bounded from below by the need to avoid 
rapid destruction and from above by the need to have eluded 
detection by direct imaging. 

We find in Section 14.71 that, if this scenario holds, the re- 
combination lines strongly overestimate the metallicities of the 
fully mixed Hii regions. The collisionally excited lines may 
also overestimate them, although in much smaller proportion 
(less than 0.1 dex in most of the cases considered in this study). 



The effect increases with ADF(0^^). In absence of any recipe 
to correct for these biases, we recommend that objects showing 
large ADFs should be discarded when probing the chemical 
evolution of galaxies. 

From our model, we do not expect a large difference be- 
tween the abundance of the fully mixed Hn region (reason- 
ably well approximated by the abundances derived from CELs) 
and the general ISM. On the other hand, ISM abundances as 
measured by absorption lines are not likely to be significantly 
affected by the presence of metal-rich droplets. On the other 
hand, the amount of O depletion of dust in the neutral ISM 
and in HII regions could be different, being more important 
in the neutral ISM than in ionized gas. Present determinations 
of the gaseous O/H in the solar vicinity from CELs in HII re- 
gions (Deharveng et al 2000) and from neutral absorption lines 
(Meyer et al. 1998, Jensen, Rachford, & Snow 2005) give sim- 
ilar values. 

In order to proceed further with this question of small-scale 
abundance inhomogeneities, which is complicated and involves 
many parameters, one needs as many observational constraints 
as possible. 

- A systematic search for ADFs for several elements (N, Ne) 
in the same regions would be welcome. If the droplet sce- 
nario is valid, ratios of ADFs for several elements in the 
same region should be in agreement with the prediction of 
current theories of nucleosynthesis in massive stars. On the 
other hand, if nitrogen is mainly produced by intermediate 
mass stars, the nitrogen abundance should not be greatly 
enhanced in the droplets we consider. The measurement of 
nitrogen recombination lines is crucial to test our scenario. 

- A systematic search for ADFs in regions of different metal- 
licities would be useful. 

- Our simulations predict that values of ADF(0'^"^) are signif- 
icantly larger than one only in specific cases. So far, ADFs 
have been measured in only a handful of Hii regions. It 
would be advisable to have a much larger number of rel- 
evant observations, and to be able to relate the values of 
the ADFs with other properties of the H ii regions and with 
the observing conditions (e.g. whether the slit encompasses 
a significant fraction of an H ii region or whether the ob- 
served fluxes are heavily weighted by zones close to the 
ionization front). It would be very useful to have an obser- 
vational census of the cases where the ADFs are close to 
one. This implies the systematic publication of intensities 
(or upper limits) of recombination lines in high signal-to- 
noise spectra. 

- Reliable measurements of ADF(O^) would also be impor- 
tant to constrain the models. At present, there are a few 
measurements of ADF(O^) in Hn regions, but they are 
highly uncertain since they rely on one or two weak lines 
from multiplet I at A ~ 111 A A (the brightest multiplet ex- 
pected to arise completely from recombination), which are 
severely affected by blends with telluric features. 

- Observational diagnostics of densities and temperatures 
from recombination lines, such as those proposed by 
Tsamis et al. (2004) or Liu et al. (2006) should be useful 
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to constrain the problem (see Garcia Rojas & Esteban 2006 
for a first attempt in this direction). 
- Another useful diagnostic is the temperature derived from 
the Balmer or Paschen jump. Recently, Guseva et al. (2006) 
have estimated the Balmer jump temperature in about 50 
metal-poor Hii galaxies. They find that the Balmer jump 
temperatures are similar to the [O iii] /14363/5007 temper- 
atures, contrary to what is found in some nearby or more 
metal-rich Hii regions (Gonzalez-Delgado et al. 2004, 
Garcia-Rojas et al. 2004). The interpretation of this result 
in the framework of the droplet scenario requires a more 
complicated modelling than that presented here, and is out 
of the scope of the present paper We just checked by hand 
on a few cases that, in our scenario, it is possible to obtain 
values of the Balmer jump temperatures equal to (or even 
slightly higher than) the [Oiii] /14363/5007 temperatures, 
depending on the fraction of [O iii] /14363 and [O iii] A5Q01 
light coming from the droplets. This question obviously re- 
quires further analysis. 

A large amount of observational work thus remains to be 
done to improve our understanding of the conditions in which 
intimate mixing occurs in the ISM. 

On the theoretical side, many questions remain unset- 
tled. For example, no firm theoretical prediction exists for the 
droplet mass spectrum, which would allow comparison with 
the constraints on droplet sizes that we derive in Section 2.2. 
Furthermore, one needs to better understand the properties of 
turbulence in H ii regions since, in order for the droplets to sur- 
vive long enough to explain the observed ADFs, it is vital that 
the regions not be too turbulent. The question of the survival of 
the droplets in the neutral/molecular ISM is another question 
that needs addressing. The photoionization models currently 
rely on several free parameters that are not well-constrained by 
observations: one needs more robust estimates of the integrated 
stellar yields as well as a better understanding of the impact of 
massive stars on the ISM. 

Given the relevance of all these issues to our understanding 
of the chemical evolution of galaxies, this subject deserves to 
be explored much further than we have been able to do in this 
paper. 
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